Information encoding in homoclinic chaotic systems 
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In this work we describe a simple method of encoding information at real time in the inter- 
spike intervals of a homoclinic chaotic system. This has been experimentally tested by means of an 
instantaneous synchronization between the laser intensity of a CO2 laser with feedback in the regime 
of Sil'nikov chaos and an external pulsed signal of very low power. The information is previously 
encoded in the temporal intervals between consecutive pulses of the external signal. The value of 
the inter-pulse intervals is varied each time a new pulse is generated. 

PACS: 05.45.Vx,05.45.Gg,05.45.-a 



Chaotic carriers are exploited to encode and transmit 
signals via two distinct approaches. The first method 
jlj consists in applying control techniques || in order to 
stabilize one of the uncountably infinite chaotic orbits 
embedded within the chaotic attractor, such that its in- 
tersections with the Poincare section can be mapped into 
a desired sequence of bits ||. The second method, that 
has been more extensively used, consists in exploiting the 
synchronization properties of coupled chaotic systems Q 
in order to modulate a chaotic carrier with a message sig- 
nal at the transmitter, and then recover the message by 
demodulating the synchronized state at the receiver. In 
the past, the latter method has been mainly focussed on 
the problem of warranting privacy in the communication 
|^| , but it has recently found applications for communi- 
cation with chaotic time-delayed optical systems ||. 

In this Letter we discuss an alternative procedure 
which makes use of pulse synchronization to encode de- 
sired messages into the interspike interval sequences of 
a homoclinic chaotic system. This way, the informa- 
tion is coded only in the time intervals at which spikes 
occur, and does not affect any geometrical property of 
the chaotic flow, thus resulting in a better performance 
against unwanted perturbations, as e.g. noise contamina- 
tion in the communication channel. We will discuss the 
performance of the proposed method in an experiment 
on a CO2 laser displaying a regime of Sil'nikov chaos. 

Sil'nikov chaos H has been observed in many systems, 
such as chemical |g] and laser || experiments. This kind 
of behavior shows striking similarities with the electri- 
cal spike trains travelling on the axons of animal neu- 
rons [jl0[ [ll]]. More generally, chemical oscillators based 
on an activator-inhibitor competition, which rule biolog- 
ical clocks controlling living rhythms, such as the heart 
pacemaker, hormone production, metabolism, etc., dis- 
play rhythmical trains of spikes with erratic repetition 
frequency as shown in Refs. (§1, ||. 

Homoclinic chaos of the Sil'nikov type normally ap- 
pears when a parameter is varied towards the homoclinic 
condition associated with a saddle focus [l2[ [l|. Its 



peculiarity consists of an astonishing regularity of the ge- 
ometric trajectory in phase space. The chaotic motion is 
characterized by the large fluctuations in the return time 
of consecutive spikes, so that only an average return pe- 
riod can be defined. Thus, an appropriate indicator of 
chaos may be the distribution of the return times to a 
given threshold, and the strength of chaos is associated 
with the amount of decorrelation between successive re- 
turns fll4|j . Such decorrelation occurs around the saddle 
focus, when the system displays a large susceptibility, i.e. 
a large response to an external perturbation Q ; we will 
exploit this susceptibility for information encoding. 

Recently, different phase locking regimes between an 
initial homoclinic chaotic signal and a sinusoidal peri- 
odic external modulation has been reported fll6|| . The 
initial homoclinic chaotic behavior can be suppressed by 
the action of a periodic external forcing, giving rise to 
different phase locking domains, (1:1, 1:2, 1:3 and 2:1), 
depending on the frequency of the applied periodic forc- 
ing. The applied forcing readjusts the evolution of the 
dynamics of the homoclinic chaotic signal in such a way 
that the return period of every orbit in phase space is the 
same. Thus, a periodic orbit is stabilized. 

Based on this synchronization, we propose an innova- 
tive and viable method of controlling instantaneously the 
temporal interval between spikes. In this way, the return 
period of each cycle changes accordingly to an external 
information. In order to get an exact control in the re- 
turn period of each orbit, it is fundamental to minimize 
the response time of the system to the forcing signal. 
For this purpose, the external perturbation is made of 
short pulses. In this way, synchronization between both 
signals is obtained almost instantaneously, thus achiev- 
ing a robust control on the individual inter-spike interval 
of a homoclinic chaotic system, and hence encoding the 
desired information in the sequence of inter-spikes. 

The experiment has been performed on a single mode 
CO2 laser with a feedback proportional to the output 
intensity (see Ref. |ll|). The control parameters, bias 
voltage and gain of the feedback, are set in a regime of 
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homoclinic chaotic behavior (see Fig. 0(a)). Figure |T](b) 
shows the phase space reconstruction of the chaotic at- 
tractor. The modulating external signal is applied, in 
an additive way, on the bias parameter with a waveform 
generator Tektronix TM5003, which is controlled by a 
real time PC board (PCI-7030/6040E) from National In- 
struments. The perturbation signal is a train of square 
waveforms. The inter-pulse intervals can be controlled, 
in real time, by means of an adequate computer program. 
Each time a new pulse is generated, this distance in time 
is changed, to encode in these time intervals the desired 
information. The inter-pulse interval is varied inside the 
domain of 1:1 phase locking for the external perturba- 
tion and the laser intensity. This time variation range 
corresponds to a frequency of repetition of pulses that 
can take values in a quite large interval. The time dura- 
tion of each pulse is 10% of the inter-pulse interval. The 
amplitude of the pulses is 4% of the value of the ampli- 
tude of the main spikes of the temporal evolution of the 
laser intensity. The laser output is recorded by a digital 
oscilloscope with a sampling time of 1 fis. 

The hrst performance index we have studied is the cor- 
relation between the different time (or frequency) inter- 
vals of consecutive pulses of the external signal and the 
corresponding return periods (or frequencies) in the tem- 
poral series of the laser intensity. This index provides in- 
formation on the amount of synchronization between the 
external forcing and the signal to control, initially in a 
homoclinic chaotic regime. Performance has been evalu- 
ated for a random uniform frequency distribution of the 
forcing signal centered in the middle value of the cor- 
responding Arnold tongue associated with the 1:1 phase 
locking, that is, centered at /o = 1.7 KHz (see Ref. [pi). 
Thus, the frequency of the external modulation take val- 
ues in the interval (/o — f s , fo + fs) provided that the 1:1 
phase locking is guaranteed, that is, for f s < 0.6 KHz. 

As shown in Fig. ||(a), the phase of the temporal series 
of the laser intensity follows almost perfectly the external 
signal, with the laser spikes adjusting in each orbit to the 
corresponding external signal. This is obtained because 
the laser parameters arc adjusted to the homoclinic chaos 
regime where the system crosses a region of large suscep- 
tibility jt5| . Figure ||(b) shows the inter-spike intervals of 
the laser intensity versus the inter-pulse intervals of the 
modulating signal. Synchronization is expressed by the 
straight line of 45 degrees. It is interesting to notice that 
synchronization does not depend on the frequency range 
of the external modulation. A better estimation of the 
error in time between the inter-spike intervals of the laser 
and the timing between pulses of the external modulating 
signal is shown in Fig. ||(c), where an histogram of these 
differences in time is plotted. Such histogram is fitted by 
a zero-mean Gaussian distribution with a standard de- 
viation of 7.8/is. In this figure we also appreciate that 
the difference between the corresponding return periods 
of the external signal and the laser intensity is always 
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FIG. 1: (a) Experimental time series of the laser intensity for 
a CO2 laser with feedback, in the regime of homoclinic chaos. 
The intensity modulation has been amplified and lifted for vi- 
sual convenience. Homoclinic chaotic behavior in a CO2 laser, 
(b) Phase space trajectory built by an embedding technique 
with appropriate delay. 



below a small value (At < 25 /is). 

From this analysis, we conclude that it is possible to 
control the laser with pulses of very low power. Thus, 
it is possible to encode in the return periods of the laser 
intensity messages of binary, ternary, quaternary nature, 
etc... The fact of encoding information in the timing be- 
tween spikes can be very useful in order to design a robust 
communication system Jl7t or to achieve a better under- 
standing in the communication phenomenon among cells 
in the central nervous system . It is well known that 
the pattern of spike times provides a large capacity for 
conveying information beyond that due to the code com- 
monly assumed by physiologists, the number of spikes 
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As an example of how to encode information in the tim- 
ing between spikes, we present the simple case in which 
a binary message is encoded in the return period of the 
spikes in the CO2 laser intensity. The homoclinic chaotic 
behavior of the laser can be used to encode a desired mes- 
sage in such a way that the bits " 1" and " 0" are easily 
identified as timing between spikes less or greater than 
a threshold value. The encoding procedure is the follow- 
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FIG. 2: (a) Experimental temporal evolution of the laser in- 
tensity driven by an external pulsed signal, (b) Inter-spike 
intervals of the laser intensity versus inter-pulse intervals of 
the forcing signal. This plot indicates synchronization be- 
tween the two variables, (c) Histogram of the time differences 
(At) between inter-pulse intervals of the external signal and 
inter-spike intervals of the laser intensity. 



ing. First of all, an external pulsed signal that modulates 
the laser intensity is synthesized with an inter-pulse fre- 
quency changed after each pulse according to a uniform 
frequency distribution centered at the mean repetition 
frequency, /o = 1.7 KHz. When the value of the fre- 
quency is larger than /o, a bit '1' is encoded in the ex- 
ternal signal; when this value is lower than /o, a bit '0' is 
encoded in the external signal. Since the error in time in 
the synchronization process between the external signal 
and the laser intensity is below 25/is, in order to avoid 
errors around the mean frequency, we consider variations 
in the frequency in the interval (/o — / s , /o — f a ) to en- 
code the bit '0', and in the interval (/o + f a , fo + fs) to 
encode the bit T'. It is obvious that for f a > 0.08 KHz 
no errors in the transmission are guaranteed, since this 
is the variation in the frequency associated with a time 
variation of 25 fis (Af — At* /q). Figure ||(a) shows the 
histogram of the inter-spike intervals of the free-running 
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FIG. 3: a) Histogram of the frequencies of the free-running 
signal before applying any perturbation, (b) Histogram when 
some information has been encoded in the signal. The infor- 
mation has been considered as a uniform distribution of bits 
T' and bits '0'. 



signal before applying any perturbation, and Figure ||(b) 
shows the same when some information has been encoded 
in the own signal. The information has been considered 
as a uniform distribution of bits T' and bits '0'. We can 
appreciate in these histograms how a left shift in time 
(right shift in frequency) is produced when the control 
is applied. This is easy to understand since the pertur- 
bation tends to eliminate the oscillations of the laser in- 
tensity around the saddle point, leading such variable to 
the rest state (no intensity). In this way, the orbit is 
shortened and the time difference between consecutive 
spikes is reduced. Notice also that the frequency range 
of the modulating signal to achieve 1:1 phase locking is 
centered around the most probable value of the frequency 
corresponding to the homoclinic behavior of the laser in- 
tensity 

The information decoding process in the laser inten- 
sity can be carried out by observing the instant when the 
laser intensity is higher than a threshold level. The fact 
that the spike amplitude is very high in comparison with 
the signal elsewhere, and the fact that the system can be 
controlled with small perturbations of very low power, 



4 



























































































































== 




























































































































i 














v 

1 



J i I i I i I i I i L_ 

5 10 15 20 25 
E b /N (dB) 



FIG. 4: BER vs. E^/No of the proposed communication sys- 
tem when additive white Gaussian noise is contaminating the 
channel. 

make it suitable for information transmission in a com- 
munication system. Since the information is contained 
entirely in the timing between spikes, channel distortions 
that affect the pulse shape will not significantly influence 
the information transmission. From this point of view, 
we analyze how robust the system is to the addition of 
a zero-mean Additive White Gaussian Noise (AWGN) 
in the communication channel. Notice that the channel 
noise is added to the transmitted signal as s(t) + n(t) in 
the transmission channel and sent to the receiver. Figure 
U shows the Bit Error Rate (BER) versus Eb/N a , where 
Eb is the energy per bit, and Nq/2 represents the power 
spectral density of the channel noise. We can see in this 
figure how the BER strongly decreases as the E^/Nq re- 
lation increases. For values of E^/No > 25 dB the BER 
is already below 10~ 3 . Own method provides a BER 
consistently lower than other recent proposals pC| [. 

In conclusion, we have shown experimentally how it is 
possible to control the inter-spike intervals of a homo- 
clinic chaotic system by applying small perturbations. 
This fact allows to encode information in the timing be- 
tween spikes. This can be useful in order to achieve a 
better understanding in the communication phenomenon 
among cells in the central nervous system as well as to 
design a robust communication system. In particular, 
we have shown a simple method to encode a binary mes- 
sage in a homoclinic chaotic signal that can be used in 
communication applications. 
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